Dynamics of bacterial flow: Emergence of spatiotemporal coherent structures 
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We propose a simple model of self-propelled particles to show that coherent structures, such as jets 
and swirls, can arise from a plausible microscopic mechanisms: (i) the elongated shape of the self- 
propelled particles with (ii) the hardcore interactions among them. We demonstrate via computer 
simulation that these coherent structures, which emerge at sufficiently high densities of particles, 
have characteristics that are similar to those observed in recent experiments in bacteria baths. 

PACS numbers: 87.18.Bb, 05.65.+b, 47.54.-fr 



o 



I 

o 
o 



> 

m 
(N 
(N 

l> 
O 

o 



X 



Living systems often exhibit complex spatiotemporal 
patterns that are characteristic of many systems driven 
out-of-equilibrium Examples range from birds flock- 
ing f^l to internal organizations inside a cell A pop- 
ulation of live microorganisms, such as bacteria like E. 
coli, suspended in an aqueous environment can provide 
a particular interesting model system to study pattern 
formation and collective motion in biology since the dy- 
namics of the individual bacteria can be directly observed 
and critical parameters such as density and activity may 
be brought under experimental control These mi- 
croorganisms continuously consume nutrients and dissi- 
pate the energy through the process of propelling them- 
selves against the frictional force exerted on them by the 
fluid ^ . With a typical size of a bacterium of the order 
of microns and a typical speed of the order of 10 /im/s, 
the Reynolds number i? ^ 1 is quite small. Yet, a large 
concentration of these microorganisms constitutes a state 
that is far from equilibrium, and exhibits self-organized 
collective motion with spatial and temporal patterns that 
are both physically fascinating and potentially of great 
biological significance [1, 0, B 0] • As a first step towards 
understanding of how suspended cells generate coherent 
motion, we identify, in this Letter, two simple but central 
ingredients - the elongated shape of the self-propelled 
particles and the hardcore interactions among them, and 
demonstrate via computer implementation of these two 
ingredients that this system exhibits coherent jets and 
swirls with characteristics strikingly similar to those ob- 
served in experiments. 

Recent experiments [1, H B Hi show that when con- 
centrated, the crowd of swimming bacteria creates arrays 
of transient jets and swirls whose size can be orders of 
magnitude larger than an individual bacterium. These 
complicated, spatially coherent structures have been ob- 
served in Bacillus subtilis colony grown on an agar plate 
0, in E. coli confined in a quasi- two-dimensional soap 
film Ml, in Bacillus subtilis at the edge of a pendent 
drop [9|. ft is estimated from direct visualization that 
these structures have a typical size of about ten times 
that of a bacterium and persist for a few seconds. More 
quantitative information about these structures may be 
extracted from a one-point passive microrheological tech- 



nique which tracks passive micron-sized beads dispersed 
in a bacterial bath of E. coli [2]. Interestingly, the mean- 
squared displacement (MSD) of these passive beads ex- 
hibit superdiffusion at short time and diffusion at long 
time. It is measured that the crossover time Tc ^ 2 s and 
the length ic = ^/{[Ax{tcW) - 10 fim provide a natu- 
ral time and length scales of these coherent structures, 
respectively. 

There have been a few simulation models and theoreti- 
cal treatments that aim to describe the above and related 
phenomena 0, [l^, [ll|, [l^] . Motivated by the patterns in 
fish schooling and bird flocking, Vicsek et al. numeri- 
cally studied a model in which each particle (modelled 
as a point) moves at a constant speed and its direction 
is determined by averaging over the directions of a large 
collection of particles in its neighborhood plus a small 
random perturbation 2]. This model in 2D exhibits a 
nonequilibrium phase transition: At sufficiently high con- 
centration, all particles spontaneously move in a single 
direction, thus breaking the rotational symmetry of the 
system. This phase transition is rationalized within a 
phenomenological dynamical xy-model [lo| . This model 
further predicts that in 2D, the MSD of a tag parti- 
cle exhibits superdiffusion at the transition. Indeed, via 
computer simulation of the bird flocking model, it was 
shown in Ref. pH that a passive bead, interacting with 
a sea of "birds" , exhibit superdiffusion as observed in 
Ref. 01 ■ However, the coherent structures appearing in 
the ordered state of the bird flocking model have a size 
scaled as the system size, whereas those observed in ex- 
periments have finite size. Note that the bird flocking 
model ignores the nematic-like ordering arising from the 
rod-like shape of a bacterium. More recently, a novel 
phenomenological macroscopic equations which take the 
ambient fluid and nematic ordering into account have 
been proposed to describe the dynamics of suspensions 
of self-propelled particles ^12\. This theory is a general- 
ization of the equilibrium hydrodynamics of nematogens 
to a non-equilibrium situation. It predicts that an or- 
dered suspension of self-propelled nematics is unstable at 
long wavelength, possibly giving rise to vortices and jets. 
Although such macroscopic equations are well suited for 
probing structures at large scales, it seems rather difficult 
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FIG. 1; Snapshots of bacterial motion at low and intermediate particle density, (a) and (b): The number of particles is 
= 2000. These two snapshots are one second apart. Each bacterium is represented as a vector whose length is the same 
as the bacterial length. There is no apparent coherent structure present. The bacteria move as if they were statistically 
independent, although some groupings occur accidentally, but they disappear quickly. The distribution is quite homogeneous 
and isotropic, (c) — 8000 — Emergence of coherent structures; note the dark stripes where a bunch of bacteria moves in 
similar directions. 



to apply them to determine microscopic mechanisms for 
the coherent structures down to scales of a few particle 
sizes. 

In this letter, we propose a new 2-D microscopic model 
that is based on two simple physical ingredients: First, 
since the ratio of the length of a bacterium to its width 
is typically ~ 5-10, we model each bacterium as an in- 
finitely thin rod with a fixed length. Secondly, we impose 
hard-core interaction among them, i.e., they cannot in- 
tersect when they are moving. Thus, our model may 
be viewed as a nonequilibrium generalization of the two- 
dimensional hard-core nematics studied in Ref. [l^. As 
we demonstrate below, our model in 2D exhibits coher- 
ent structures which emerge at sufficiently high densities, 
even though the bacteria in our model have no "will" , in 
the sense that they do not seek out at what directions 
their neighbors are swimming, they do not communi- 
cate with each other via chemical signaling, nor do they 
interact hydrodynamically. Rather, the simple physical 
picture emerging from this study is that, with the con- 
straints of the elongated shape of the needle and the ex- 
cluded volume interaction, a self-propelled bacterium can 
find a new location only where it can fit into. This forces 
the needles to align locally and swim together in the form 
of swirls and jets. Our simulation shows that (i) typical 
sizes of these coherent structures are consistent with ex- 
perimental observations, and (ii) the density fluctuations 
are characterized by non-Poisson statistics. 

We now detail our model of bacterial bath and de- 
scribe briefly our simulation method. We first deposit 
N needles of length £, whose center-of-mass coordinates 
Tj — ixi,yi), i = 1,...,N, are drawn at random from 
a uniform distribution in a 2D box with the length L. 
The orientation of each needle can be characterized 
by the angle with respect to the y-axis. This angle is 
chosen at random in the range of [0, 27r] from a uniform 
distribution. At this stage, there may be many bacteria 



intersecting with each other. To impose the hard-core in- 
teractions, we must first obtain a non-intersecting state. 
This is accomplished by moving each needles by a certain 
distance and rotating its direction until no needles are 
intersecting. This state constitutes an acceptable initial 
state. We use the criterion in Ref. [l^ to test whether 
two needles are intersecting: Two needles i and j inter- 
sect if and only if both gi and gj are negative, where the 
quantities gi and gj are defined as 

5. ^ (r,,.v,)2-(^/2)2[l-(v,.v,)^], (1) 
g, ^ (r,, .v,)2-(^/2)2[l-(v,.v,f], (2) 

with rij — r-i — Vj and and Vj being unit vectors per- 
pendicular to hi and iij, respectively. We employ peri- 
odic boundary condition in both directions. The bacte- 
rial motion is simulated with the following dynamics: a 
particle is chosen to move to a new position while remain- 
ing particles are fixed. Note that this asynchronous up- 
dating avoids complication associated with many-body 
collisions, which would be difficult to resolve computa- 
tionally for the hardcore interactions. Thus, at each time 
step, an individual bacterium moves by a fixed distance d 
and simultaneously its direction changes by a small angle 
6, which is chosen uniformly at random from the inter- 
val [—77/2,77/2]. If the chosen bacterium intersects with 
another one, in order to satisfy the hard-core constraint, 
we immediately discard the move and start a new trial 
move: we pull it back to the previous (original) posi- 
tion, and change its direction by another random angle 
6. If the number of such trials exceeds a preset number 
A^triai ~ 40 14], then we increase 77 by a small amount to 
a new rj' = a-q, (say, a — 1.02) and repeat the trial until 
this particular particle can make its move. Then we reset 
rj to its default value and update the next particle. We 
emphasize that there is no preferential searching involved 
in a sense that the bacterium is merely trying to find a 
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FIG. 2: Snapshots of bacterial motion at the high particle density (A'^ = 12000). See Fig. ^ for the graphic representation 
convention. Dense needle regions form coherent structures of jets and swirls. Some of these structures can be seen to persist 
for seconds as shown here. (The time between each snapshot here is 1 s.) Note the temporal evolution of the shape and size of 
these structures. 




FIG. 3: Snapshots of velocity fields, A'^ = 12000, (see text). Long jets and large swirls are persistent. The size of jets is 
comparable to the box size; swirls can have long life-time (around 8—10 seconds). The time between each snapshot here is 1 
second with a total duration of two seconds shown.) 



new location where it can fit into. This procedure to up- 
date the position and direction of the needle is repeated 
for all particles sequentially at each time step. We carry 
out our simulations, and discard the first ~3000 moves 
as transient. For all the simulations reported below, we 
have chosen i = 0.04, L = 1, d = 0.001, and N ranges 
from ~ 10^ to ~ lO''. The parameter 77 determines how 
widely a bacterium can change its direction, and its de- 
fault value is 77 = 0.01. In order to compare with physical 
values, each time step corresponds to ^ 0.01 s. Roughly, 
this equivalently corresponds to the situation where a 
bacterium of length of a few microns swims at a typical 
speed of a few fim/s 0]. 

We now turn to our main simulation results. Depicted 
in Figs. [TJ and [2] are snapshots of the needle configura- 
tion for N = 2000, 8000, and 12000, respectively. As 
evident in these figures, different densities of particles re- 
sult in different structures in their collective motion. For 
low densities, around N ~ 10^ there are no spatiotempo- 
rally coherent structures, as shown in Fig. [TJa) and (b). 
Bacteria rarely collide. However, as the concentration 



of the particles reaches the values of 4000, some 

coherent structures start to emerge with a short persis- 
tent time {t ^ Is). At = 8000, these coherent struc- 
tures become more prominent and one can easily identify 
long jets and small swirls (Fig. mc)). These structures 
can persist for a few seconds. For similar densities, the 
same structures are observed experimentally 0, . As 
the particle concentration increases, the typical size of 
the coherent structures and the frequency of their occur- 
rence also increase. The highest density we consider is 
A^ = 12000. In this case, jets of the size of the box appear 
and almost all particles are involved in collective motion 
[Fig. [2]. The evolution of the velocity field is shown in 
Fig. [31 where the velocity field is the average vector by 
summing bacterial vectors locally with the averaging area 
being £'^. Figure [3] clearly demonstrates the persistence of 
coherent structures for seconds in the velocity field. At 
a fixed particle concentration, the size of swirls clearly 
depends on rj. As rj increases, the typical size of swirls 
decreases. This is because particles in the swirls at each 
time step explore wider angles and are able to escape 
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FIG. 4: Particle density fluctuation characterized by Cy . 
Tcv = 0-1 s. For a fixed A'^, we plot here the value of Cy^ 
for each independent run. For example, at A'^ = 8000, there 
were three runs, thus three different values of Cy^ . (Instead 
of averaging over these runs, we display this way to indicate 
the statistical spread.) 



from the region of the swirl. For r/ ~ 0.01, the size of 
swirls is about 5-7 times the lengths of the bacterium at 
N ~ 4000-8000. One of the temporal characteristics of 
jets or swirls is their persistence time. Experimental ob- 
servations indicate that the persistence time for swirls 
is on the order of several seconds, while jets can persist 
for longer time. Our simulations yield similar results for 
swirls and jets. Incidentally, for initial conditions with 
all parallel needles, our simulations show that this initial 
long-ranged order is quickly destroyed but the coherent 
structures of jets and swirls again develop as the sys- 
tem evolves. We emphasize that the hard-core interac- 
tion among elongated self-propelled objects underlies all 
these structures in the simulation. The physical picture 
emerges from this model is that the excluded-volume in- 
teractions effectively reduces t], which forces particles to 
align locally. This mechanism appears to be consistent 
with experimental observations in that bacteria tumble 
(hence, changing directions) less when concentrated 0. 

To further quantify the transition from the homoge- 
neous, isotropic flow to highly organized coherent mo- 
tions of bacteria, we measure the coefficient of variance, 
Cy , from our simulation to describe the coherence of the 



motion, where Cy = a'^/{n), (n) is the average number 
of bacteria appearing in a small square with an area 1/16 
over a fixed time interval Tcy , and cr„ is the number fluc- 
tuation over Tcy ■ Figure 3] displays that Cy ^ is nearly 
one at the low density limit (Note that Cy = 1 for a 
Poisson statistics.). Then, Cy^ has a rapid rise around 
N ~ 10"*, signifying a drastic reduction of fluctuations 
in the particle density. Hence, the onset of the coherent 
motion of the bacterial flow. 

In conclusion, we have proposed a possible mechanism 
for the emergence of coherent structures in a bacterial 
flow: they arise in 2D from their elongated shapes with 
the hard-core interaction. We have shown that the model 
indeed gives rise to persistence of jets and swirls at suf- 
ficiently high density of needles, with comparable spa- 
tiotemporal scales as observed in experiments for bacte- 
ria at similar densities. When the bacterial concentration 
becomes sufficiently high, one may suspect that hydro- 
dynamic interactions of solution become important. As 
noted before, our simple model does not taken into ac- 
count hydrodynamic interactions. In general, inclusion of 
hydrodynamic interactions for elongated objects is com- 
putationally involved. However, we have performed a 
simulation of spherical particle motion following method 
in Ref. inside a 3D unit box, taking into account of 
the lubrication effect only. In this setup, we did not ob- 
serve any jets or swirls, which means that hydrodynamics 
alone probably cannot organize spherical particles into 
large-scale coherent motions. This is indicative of the 
fact that the elongated shape together with hardcore in- 
teraction in 2D plays a primary role for the emergence of 
coherent structures even in the presence of hydrodynam- 
ics. We believe that at high densities the formation of 
the coherent structures due to our mechanism is robust 
qualitatively, with possible quantitative modifications by 
details of hydrodynamic interactions. 
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